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I.  INTRODUCTION 


This  is  the  fourth  quarterly  technical  report  for  contract 
No.  F30602- 72-C-0305 .  The  object  of  the  program  is  to  provide  theoret¬ 
ical  studies  of  selected  topics  associated  with  light-beam  propagation 
through  a  turbulent  atmosphere.  The  present  problem  of  interest  has  to 
do  with  atmospheric  imaging  and  image  restoration.  It  is  ar,  examina¬ 
tion  of  the  effectiveness  of  point  and  other  types  of  reference  sources 
for  image  restoration  in  the  cases  where  the  restoration  reference 
source  and  object  are  not  necessarily  at  the  same  range  and  positior  a* 
the  object.  The  general  area  is  of  particular  interest  to  those  in¬ 
terested  in  taking  precision  photographs  of  earth  satellites  with 
ground-based  telescopes.  It  may  also  be  of  interest  to  those  consider¬ 
ing  doing  long  distance  imaging  of  terrestrial  objects. 

The  work  presented  in  this  report  is  based  on  a  familiar  procedure 
for  image  restoration:  the  Fourier  transform  procedure  (Harris,  1066 ) 
There,  one  images  simultaneously  a  turbulence  degraded  object  and  a 
turbulence  degraded  point  reference  source.  If  every  ooint  in  the 
object  and  the  point  reference  source  are  degraded  exactly  the  same, 
then  the  quotient  of  the  degraded  object  spatia1  spectrum  and  the 
reference  point  spatial  spectrum  is  the  restored  image  spatial  spectrum. 
This  is  the  ideal  case.  The  investigation  reported  herein  considers 
he  situation  where  the  reference  source  is  differert  in  location  and 
has  slightly  different  degradation  from  that  of  the  ideal  situation. 

The  situation  of  interest  is  shown  in  Fig.  1.  We  consider  two 
point  sources ,  possibly  at  different  distance1"  from  an  imaging  telescope 
The  light  from  the  two  sources  travels  through  a  turbulent  atmosphere 
and  is  imaged  on  two  image  planes  appropriate  to  their  distances.  One 
then  records  the  instantaneous  intensity  patterns  and  uses  them  with 

the  appreciate  transverse  size  scaling  in  the  Fourier  transform  restora 
tion  procedure. 

The  restoration  procedure  in  this  situation  is  expected  to  have 
limitations.  The  situation  becomes  complicated  if  1)  the  object  points 
and  reference'  point  are  not  at  the  same  transverse  location  so  that 
rays  from  the  two  partly  intersect  areas  of  different  index  fluctua¬ 
tions,  and  2)  the  object  and  point  reference  sources  are  at  different 
distances  from  the  receiver  so  that  refractive  index  fluctuations  of  a 
given  size  experience  different  magnifications  by  the  time  they  have 
arrived  at  the  receiver.  J 

.  ..  ?"e  way  t0  overcome  the  limitation  arising  when  the  object  is  at 
infinity  and  the  reference  source  is  at  a  finite  distance  from  the 
receiver  is  to  use  a  collimated  rather  than  a  power  reference  source. 

This  would  duplicate  the  point  source  at  infinite  distance  and  should 

work  if  the  degrading  turbulence  is  between  the  reference  source  and 
the  receiver. 
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The  problem  considered  in  these  calculations  is  to  assess  the 
effects  of  displacer  ant  ol  the  reference  source  and  to  determine  the 
feasibility  of  using  the  plane-wave  source  as  a  possible  remedy  in 
overcoming  the  effects  of  longitudinal  distance  difference. 


A  model  is  chosen  which  allows  the  problem  to  be  g-eatly  simpli¬ 
fied.  The  model  is  characterized  by  the  assumption  that  the  effect 
of  the  atmospheric  turbulence  on  the  wavefront  of  the  light  arriving 
at  the  imaging  telescope  aperture  can  be  represented  by  deviations  from 
the  non-turbulent  situation.  These  deviations  are  described  nathe- 
matically  in  terms  of  a  series  of  polynomials  in  transverse  positior, 
r,  which  indexes  the  deviation  at  any  given  aperture  position.  Thus 
we  reprerent  the  phase,  $,  of  the  light  arriving  at  position  r  in  the 
aperture  from  a  given  source  in  the  form 


*(r)  =^o  (r)  +  <M") 


where 


*i(r)  =  [  an  Fn(r) 
n 


and  where  4>0  represents  the  phase  in  the  ideal  (non-turbulent)  case  and 
<i(r)  represents  the  atmospheri cally  induced  fluctuations.  The  poly¬ 
nomial  series  similar  to  that  used  by  Fried  (Fried,  1967)  will  be  given 
later. 


The  approach  to  the  problem  has  two  parts,  both  of  which  can  be 
stated  in  terms  of  the  quantities  in  Eqs.  (1).  In  the  first  part,  the 
Fourier  transform  restoration  procedure  is  carried  out  using  a  parti¬ 
cular  model  for  the  object  and  reference  input  plane  field.  In  the 
model  these  fields  are  represented  using  one  or  more  terms  in  the 
series  in  Eq.  (lb),  thus  representing  particular  wavefront  deviatirns. 
Indeed  for  the  case  considered,  one  polynomial  (the  same  one)  is  used 
for  both  the  object  and  reference  beam  phases.  The  two  wavefronts 
differ  in  the  size  of  the  coefficients  an  and  an'  multiplying  the  poly¬ 
nomials;  i.e.,  they  have  the  same-shaped  wavefront  deviation  but  one 
is  more  bent  than  the  other.  The  quality  of  the  restoratinn  is  then 
investigated  as  a  function  of  the  difference  between  the  coefficients. 
This  allows  a  limit  to  be  set  on  the  size  of  the  coefficient  difference 
before  the  restoration  procedure  breaks  down. 


The  coefficient  difference  is  further  related  to  the  atmospheric 
parameters.  Specifically,  the  root-mean-square  coefficient  difference 
is  calculated  as  a  function  of  object-to- reference  source  separation, 
range,  wavelength,  turbulence  strength,  and  receiving  plane  aperture 
size,  etc. 


The  combination  of  the  image  processing  and  atmospheric  calcula¬ 
tions  for  the  coefficient  difference  then  allows  one  to  calculate  an 
index  of  successful  restoration  for  a  given  set  of  atmospheric  parameters 
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It  is  assumed  that  the  signal-to-noise  ratio  is  sufficiently  large 
so  as  to  allow  complete  recording  of  all  pertinent  images. 

In  a  previous  report  (  3432-4)  the  analysis  of  the  restored  image 
spectrum  in  the  first  part  of  the  analysis  was  started.  The  work 
presented  in  this  ,eport  extends  the  part  of  the  effort  concerned  with 
image  restoration  studies.  The  work  presented  also  represents  the 
beginning  of  the  other  part  of  the  effort,  the  atmospheric  part. 

The  material  presented  in  this  report  is  by  no  means  complete  but 
summarizes  work  in  progress.  It  draws  on  two  other  reports:  3432-5 
which  covers  the  restoration  studies  in  more  detail,  and  3432-7  which 
presents  the  details  of  the  atmospheric  calculation.  Further,  neither 
of  these  reports  represents  a  completed  work,  but  both  present  detailed 
calculations  from  which  preliminary  results  may  be  drawn. 

In  section  II  of  this  report,  the  Fourier  transform  restoration 
procedure  is  applied  to  images  formed  with  wavefront  image  degradation 
determined  by  one  or  another  of  the  wavefront  polynomials.  An  expres¬ 
sion  for  the  restored  image  spectrum  is  stated  under  the  assumption 
that  only  the  first  six  polynomials  in  the  series  have  significant 
effect.  In  section  III  three  criteria  for  restored  image  evaluation 
are  presented.  Section  IV  contains  examinations  of  restored  images 
using  the  three  image  criteria  as  the  difference  in  polynomial  coeffi¬ 
cients  is  increased.  This  leads  to  possible  estimates  of  upper  limits 
on  the  wavefront  differences.  Section  V  contains  an  analysis  of  the 
wavefront  deviation  determined  by  particular  atmospheric  conditions. 

The  effects  of  standard  ground  layer  turbulence  are  considered,  as  are 
those  of  a  strong  turbulent  layer.  The  summary  is  in  section  VI. 

II.  IMAGE  PROCESSING  EFFECTS 

In  this  section  we  examine  the  results  of  applying  the  Fourier 
transform  restoration  procedure  to  images  degraded  by  the  addition  of 
the  particular  polynomial  shape  to  the  incoming  ideal  electric  field 
phase.  The  restored  images  are  compared  using  several  different  cri¬ 
teria  to  give  information  about  allowable  wavefront  differences. 

The  approach  follows  that  introduced  in  the  previous  quarterly 
r  port  (3432-4).  The  electric  fields  arriving  at  the  input  aperture 
from  both  the  object  and  reference  point  sources  are  assumed  to  have  the 
form 


(2)  E  =  E0  exp(j<j>(r)) 


where  <*>(7)  is  given  by  Eqs.  (1).  The  orthonormal  polynomials  used  are, 
up  to  the  second  order, 
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(3a) 

(3b) 

(3c) 

(3d) 

(3e) 

(3f) 


r  _  /3  X 

F  ,,ijL 

3  a  a/2 

F4  -  J-fjf  (x2  +  y2  -  (a2/6))/(a/2)2 

Fs  -  j  |  J  |  (x2-y2)/(a/2)2 

F6  =  (3/ a)  (xy)/(a/2)2  . 


These  polynomials  are  chosen  so  as  to  be  orthonormal  over  a  square 
aperture  of  side  "a".  It  is  assumed  that  the  major  effect  of  the 
atmosphere  is  represented  by  the  first  six  polynomials,  an  assumption 
born  out  (at  least  in  the  mean  square  sense)  by  out  calculations  and 
by  tried. 

The  first  step  in  using  the  polynomial  procedure  is  to  calculate 
an  expression  for  the  restored  image  spectrum  assuming  the  input  light 
field  described  in  Eq.  (2).  This  involves  first  finding  the  field 
correlation  functions  of  the  object  and  reference  fields.  These  repre¬ 
sent  the  spatial  spectra  of  the  associated  object  and  reference  image 
intensities.  The  quotient  of  i-hese  two  spectra  is  the  spatial  spectrum 
of  the  restored  image.  The  expression  for  the  restored  image  spectrum 
is  then  (3432-4), 


(4a)  W0(k) 


w '  (~) 


^xxr  x+^xy^/  ^yy'y, +^x/xj 


sin[((ka/d)-|Kx|)(BxxKx*BxyKv)j  sin|((ka/d)-|xy|)(ByyCy+BxyKx)' 
sfnj(  (ka/d)-|<xj){BxxKx+BxyKy)|  sin  [( (ka/d)-|xy|  XB^y-B^)  | 


ka 

d 


<  K  < 

X 


+ 


d 


ka 

d 
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where  the  quantities  B  ,  Bxv,  Bvv  and  4  are 
deformation  parameters,  a2  through  a6  by0 


related  to  the  wavefront 


(4b)  *0  =  2/J(d/a2k)ta2<x  +  a^y) 

(4c)  Bxx  =  6^72  (a4  +  a5)  ( d/2k ) 2/ ( a/2 ) 3 

(4d)  Bxy  =  Byx  =  6a6(d/2k)2/(a/2)3 

(4e )  Byy  =  6^5/2  (a4  -  a5)  ( d/2k ) 2/ ( a/2 ) 3  . 


The  p 'i mes  indi cate  the  coefficients  for  the  reference  beam,  the  un- 
primr  d  vari ables  are  for  the  object  beam.  Of  the  other  quantities, 
k  =  2tt/a,  where  A  is  the  light  wavelength,  "a"  is  the  size  of  the  square 
telescope  receiving  aperture,  "d"  is  the  telescope  image  distance,  and 
kx  and  Ky  are  spatial  frequency  components. 


iS  no^e<^  That  if  the  object  and  reference  wavefronts  are 

identical;  i.e.  an  =  an  1 ,  then  all  the  primed  and  unprimed  quantities 
are  equal  respectively,  and  W0(<)  is  equal  to  unity. 


It  is  also  noted  that  the  linear  wavefront  tilt  terms  a?  and  a? 
produce  only  translation  of  the  restored  image  since  they  only  appear 
in  the  spatial  spectrum  as  the  factor  eJ(*o-$o  with  an  extended 
object  these  would  move  points  around,  but  not  degrade  resolution  For 
this  reason  effort  is  concentrated  on  estimating  the  effect  of  the  qua¬ 
dratic  wavefront  distortion  terms  indexed  by  the  B.  .  and  B-1-. 

1  J  1  J 


The  expression  in  Eq.  (4)  is  evaluated  using  selected  values  for 
the  parameters.  It  is  digitally  Fourier  transformed  and  the  resultinq 
two-dimensional  images  are  examined  using  several  different  criteria. 


III.  IMAGE  EVALUATION  CRITERIA 


Three  different  criteria  were  used  to  evaluate  the  restored  images 
They  are  the  restored  image  Strehl  ratio,  the  normalized  integral  scale’ 
widtn,  and  the  mean  square  difference  from  the  ideal  restored  image. 

These  criteria  give  respectively  measures  of  the  restored  image  contrast, 
the  restored  image  resolution,  and  the  intpnsity  difference  from  the 

m^n  square  error  is  deluded  because  it  is  conceivable 
tnat  the  Strehl  ratio  which  is  proportional  to  the  integral  over  spatial 
frequency  space  of  the  restored  image  spectrum  might  have  a  reasonable 
value  even  though  the  restored  image  spectrum  was  fluctuating  wildly 
but  having  nearly  equal  positive  and  negative  parts.  This  would  give  a 
reasonable  Strehl  ratio  but  a  poor  restored  image.  Similar  behavior  could 
occur  for  the  integral  scale. 
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Numerical  values  were  chosen  to  indicate  whether  an  image  was 
properly  restored  according  to  these  criteria.  Specifically,  for  the 
Strehl  ratic  and  the  normalized  integral  scale  width,  they  were  "three 
db"  points.  Thus  if  the  Strehl  ratio  and  normalized  integral  scale  of 
the  restored  image  were  between  0.5  and  2.0,  the  image  was  termed 
acceptably  restored.  Similarly,  a  mean  square  deviation  of  0.1  or  less 
was  termed  acceptable. 

IV.  IMAGE  PROCESSING  STUDIES  RESULTS 

The  details  of  the  calculation  of  the  restored  images  and  the 
evaluation  of  the  three  measures  of  the  image  quality  are  given  in 
another  report  (3452-5).  In  order  to  understand  the  detailed  behavior 
of  these  three  quantities,  it  was  necessary  to  examine  in  detail  the 
locations  of  the  poles  and  zeros  of  the  restored  image  spectrum. 

Indeed,  the  one  most  significant  factor  in  determining  the  behavior  of 
these  quantities  was  the  change  of  these  poles  and  zeros  with  the  size 
of  the  a^,  the  wavefront  polynomial  coefficients.  These  factors  are 
covered  in  detail  in  report  3432-5. 

T .«?  ^sults  o*  the  restored  image  evaluation  are  presented  in  a 
particular  form.  Implicit  is  the  assumption  that  both  the  object  and 
reference  fields  are  assumed  to  be  perturbed  by  the  sa'ie  shape  functions 
They  differ  in  the  amount  of  one  function.  The  wavefront  deformation  is 
represented  by  the  coefficients  an,  and  an'  for  the  object  and  reference 
values,  respectively.  A  particular  value  is  chosen  for  one  of  the  co¬ 
efficients,  say  a^  and  a  slightly  different  value  is  chosen  for  am' , 
the  difference  being  designated  ny  Aam  where  Aam  =  aJ  -  am.  The  image 
evaluation  criteria  are  then  plotted  as  a  function  of  Aam.  Thus,  in  the 
case  where  Aam  =  0,  there  may  be  considerable  degradation,  i.e.,  large 
values  of  the  an,  but  still  perfect  restoration. 

One  case  is  shown  here,  a  case  where  all  but  a4  aid  ai/  are  zero, 
these  coefficients  representing  focusing  errors.  The  case  chosen  has 
an  object  wa/efront  deviation  of  three-quarters  of  a  wavelength  at  the 
aperture  edge  on  axis,  a  value  sufficiently  large  so  as  to  give  observ¬ 
able  degradation. 

The  restored  image  Strehl  ratio,  integral  scale  and  mean  square 
errors  are  shown  in  Figs.  2a,  2b  and  2c,  respectively,  plotted  as  a 
function  of  wavefront  difference  at  the  aperture  edge  on  axis.  There 
are  several  items  of  interest  in  these  figures.  First  consider  the 
Strehl  ratio  as  shown  in  Fig.  2a.  We  see  that  when  there  is  no  wave- 
front  difference,  Aa4  =  0,  the  Strehl  ratio,  Sr,  equals  one  as  it  should 
For  values  of  Aai*  between  -0.003  and  +0.005  there  is  acceptable  restor¬ 
ation.  For  values  of  Sr  outside  that  region  there  are  alternating 
regions  of  unacceptable  restoration,  indicated  by  the  cross-hatched 
areas,  and  acceptable  restoration.  It  is  worthwhile  pointing  out  the 
fact  that  the  image  quality  does  not  slowly  deteriorate  and  then  stay 
bad,  but  has  alternating  bad  and  good  regions.  Similar  behavior  occurs 
for  the  integral  scale  and  mean  square  error. 
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restored  image  Strenl  ratio  versus  wavefront  difference  for  object-wavefront 
of  three-quarters  of  a  wavelength. 


mags  integral  scale  width  versus 
of  three  quarters  of  a  wavelena 


Also  of  interest  is  that  the  restored  image  spectrum  and  the 
restoration  process  depend  on  the  number  of  wavelengths  deviation  of  the 
image  wave  itself.  That  variation  also  is  not  monotonic  in  nature. 
Further,  graphs  given  in  report  3432-5  show  that  for  0.5  wave¬ 
lengths  deviation  at  the  aperture  edge  on  axis,  there  is  difficulty  in 
restoring  the  image  because  of  the  number  of  poles  and  the  way  they 
move  with  respect  to  the  sampling  points.  The  case  is  worse  than  for 
the  one  and  one-half  wavelengths  case  shown  here.  On  the  othir  hand, 
the  case  where  there  is  a  half-wave  variation  for  the  image  plane 
wavefront  on  axis  is  relatively  easy  to  restore. 

It  should  be  noted  that  the  data  presented  in  Figs.  2  was  obtained 
by  a  dogmatic  application  of  the  Fourier  transform  procedure.  No 
attempt  was  made  to  smooth  or  clip  the  restored  image  spectrum.  There 
is  a  significant  question  as  to  whether  the  discrete  representation  used 
has  a  sufficiently  fine  point  spacing  to  faithfully  represent  the  func¬ 
tion  without  some  smoothing,  and  meaningfully  with  smoothing. 

To  summarize,  the  Fourier  transform  restoration  procedure  has  been 
applied  to  simple  point  object  and  reference  images  with  simple  varia¬ 
tion  of  the  wavefront  from  ideal.  Using  the  Fourier  transform  procedure 
with  no  restored  image  spectrum  smoothing,  an  examination  was  conducted 
to  see  how  much  difference  could  be  allowed  between  the  object  and 
reference  wavefronts  before  restoration  became  impossible.  Typical 
nunfcers  are  .003  wavelengths  difference  at  the  edge  to  assure  restor¬ 
ation,  and  larger  than  0.1  wavelengths  if  one  accepts  alternating 
regions  of  restorable  and  non-restorable  images. 

V.  ATMOSPHERIC  EFFECTS 

In  this  section  a  study  is  made  of  the  size  of  the  maximum  wave- 
front  deviation  of  the  spherical  and  other  components  due  to 
atmospheric  refractive  index  effects.  The  input  information  comprises 
physical  parameters  such  as  ranges,  separation  of  the  point  sources, 
and  variation  of  the  refractive  index  fluctuations.  The  questions 
asked  relate  to  the  root-mean-square  values  of  the  atmospherically 
induced  variation  of  the  coefficients  in  the  wavefront  polynomial  series 
for  three  source  configurations  and  for  ground  level  turbulence  and  high 
altitude  turbulent  layers.  The  calculation  uses  the  (reasonable) 
assumption  that  phase  effects  are  predominant,  and  neglects  amplitude 
effects.  Results  will  be  summarized  here,  and  the  details  presented  in 
a  different  report  (3432-7). 

The  physical  situation  considered  in  the  analysis  is  shown  In 
Fiq.  3.  There  we  see  two  point  sources  located  at  points  r  and  r+p. 

The  phase  variations  from  ideal  of  the  light  arriving  at  the  input 
aperture  from  the  two  sources  are,  respectively,  $,(r,,r)  and  i,(ri,r+7) 
given  by  1 


. * 


. . . . 


(6a)  4-i(r,rj)  =  J  an(rj)  Fn(r) 

and 

(6b)  +i(r,r1+(7)  =  £  Sp’^+p)  Fn(r) 

The  light  from  the  two  sources  may  travel  through  regions  of  different 
atmospheric  index  fluctuations,  and  therefore  have  different  series 
represented  by  the  an(rj)  and  ap'^+iD. 

The  quantity  of  interest  is  the  root-mean-square  coefficient 
di  fference. 


(7)  *Da(F)  =  /<|an'(F1+J)  -  a„(r1) |2) 
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where  the  brackets  represent  ensemble  averaging. 

The  input  formulae  used  for  the  calculation  of  4)a(p)  are  the 
definition  of  the  an  based  on  the  fact  that  the  polynomials  Fn(r)  are 
orthonormal  over  a  square  aperture  of  side  "a" 


(8)  an  =  j|Fn(r)  (F)  W(r)  dr 

(W(r)  is  the  aperture  function  equal  to  1  inside  the  aperture  andO  outside 
it),  the  expression  for  the  phase 

L 

(9)  =  k  (  n(r')  dz 

0 

and  the  expression  for  the  index  structure  function. 


(10) 


Dn  (r'-r")  =  Cn2  Ir'-r"!2^  . 


Combining  Eqs.  (7),  (8),  (9),  and  (10)  gives  the  basic  expression  for 

_  I  )  % 

(11a)  ^Da(p)  =  |  jjdr  dr*  Fn(r)  Fn(r')  W(r)  W(r')  D^r-r'.p)  j> 


where 
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(lib)  D^r-r'.p)  =  i  •'  d2j  dz2L((p1-p2,)f  Zl-Zp) 

0  0 


+  Dn^P2_Pl  ’  *2l-22  )  "  C,n(pl"p2’2l"z2)  "  Dn^p2-P2  1  »2l-22)  * 


The  integrals  in  Eq.  (11)  are  simplified  and  then  evaluated  using 
a  digital  computer.  The  details  are  presented  in  report  3432-7;  only 
typical  results  w~ll  be  presented  here.  Such  results  are  shown  in 
Fig.  4  where  we  see  /Da4(p)  graphed  as  a  function  of  p  for  point  objects 
at  20  km  and  150  miles  and  point  reference  at  20  km,  and  for  two  turbulent 
conditions.  The  two  bottom  curves  show  the  effects  of  surface  level 
turbulence  (Cn2(z)  =  (z/Zq)"^3)  for  the  two  different  objects  and  the  two 
top  curves  show  the  effects  of  a  high  altitude  turbulence  layer  at  17  km. 
These  curves  apply  specifically  for  coefficient  a.  ;  however,  similar 
curves  are  shown  for  a5  and  a6  in  report  3432-7. 

The  strength  of  the  turbulent  iayer  was  chosen  so  as  to  provide  the 
same  value  for  a  plane  wave  wave  structure  function  at  the  ground  for  a 
point  source  at  infinite  range  and  turbulent  layer,  as  due  to  ground 
level  turbulence,  and  point  source  at  infinity.  The  calculations  were 
made  forx=  0.5  y  and  Cn2  =  2.1  x  10“16  at  one  meter  height. 

The  curves  exhibit  some  expected  behavior.  For  the  cases  where 
point  object  and  reference  are  both  at  20  km  height,  Ah (P)  goes  to  zero 
as  p  becomes  small.  This  is  expected  since  light  from  both  (point) 
sources  ends  up  traversing  regions  of  the  same  refractive  index  fluc¬ 
tuations.  For  the  case  where  one  sourr*  is  at  a  different  height  than 
the  other,  there  is  a  residual  value  f  /Daif(0)  even  when  p  =  0.  This 
is  because  a  given  index  fluctuation  magnifies  liqht  differently  from 
sources  at  different  distances,  and  there  is  complete  cancellation  only 
for  very  low  level  turbulent  fluctuations.  It  is  also  noted  that  the 
effect  of  the  high-level  turbulent  layer  is  significantly  greater  than 
that  due  to  the  ground  level  turbulence.  That  is,  a  small  separation 
of  rays  from  the  two  sources  produces  a  much  greater  effect  in  inter¬ 
acting  with  high  level  index  fluctuations  than  with  ground  level 
fluctuations,  due  to  the  wavefront  magnification  of  spherical  waves. 

The  curves  can  be  useful  in  deciding  whether  restoration  is  possible 
in  a  given  situation.  For  example,  if  a  value  of  0.01  wavelengths  dif¬ 
ference  can  be  tolerated  for  the  wavefronts  then  a  verv  close  spacinq 
between  reference  and  image  must  be  obtained  or  there  must  be  weaker 
turbulence.  If  a  value  of  0.10  is  acceptable,  then,  with  the  value  of 
cn  chosen,  a  source  separation  of  one  aperture  size  is  allowed. 
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SEPARATION*  RH0 

Graph  of  RMS  wavefront  difference  versus  object-reference  source  distance. 


One  other  point  is  indicated.  That  concerns  the  relative  ease  of 
restoration  for  objects  at  the  same  height  as  the  reference,  as  for  an 
object  much  higher  than  the  reference.  Fig.  3  indicates  that  for  object 
and  reference  separated  by  more  than  one  aperture  size,  the  ease  of 
restoration  is  almost  the  same,  the  values  of  vCjrT  lying  within  a 
factor  of  two  of  each  other  for  p/a  >  i. 

To  summarize  this  section,  the  root-mean-square  coefficient  differ¬ 
ence  induced  by  turbulent  index  fluctuations  has  been  calculated. 

Results  indicate  that  high  altitude  turbulence  effects  are  somewhat 
stronger  than  those  due  to  ground  layer  turbulence.  Further,  the  ability 
to  restore  is  about  the  same  for  an  object  at  the  reference  source 
heigh,  as  it  is  for  an  object  much  higher  than  the  reference  source  for 
very  reasonable  separations  between  object  and  reference.  RMS  coefficient 
values  between  10"  and  10"1  wavelenqths  are  to  be  expected. 

VI .  SUMMARY 

Summarizing  the  results  of  the  report,  the  problem  of  using  the 
Fourier  transform  restoration  nnocedure  to  restore  atmospherically 
degraded  images  was  divided  -*iito  two  parts,  a  study  of  the  restoration 
procedure,  and  a  study  of  atmospheric  effects.  The  two  parts  were 
joined  by  choosing  a  model  which  assumes  the  phase  front  of  an  atmos¬ 
pherically  degraded  light  wave  can  be  represented  as  a  truncated  poly¬ 
nomial  series.  The  effects  of  a  wave'ront  shape  represented  by  one 
term  of  the  series  was  examined,  in  this  case  the  shape  representing 
atmospherically  defocused  images.  Preliminary  results  indicate  that  the 
e.tects  of  wavefront  differences  between  object  and  reference  beams  are 
not  simply  described,  being  an  alternating  function  of  wavefront  differ¬ 
ence  and  not  a  monotonic  function.  Wavefront  differences  of 
ess  than  .003  wavelengths  assure  restoration  effectiveness,  while  much 
larger  wavefront  differences  may  allow  restoration. 

As  far  as  atmospheric  effects  are  concerned,  root-mean-square 
wavefront  differences  of  .001  to  .01  wavelengths  are  reasonable  for 
nedium  strength  ground-level  turbulence.  High  altitude  turbulence 
layers  and  strong  ground  level  turbulence  may  make  restoration  impossible 
A  p°int  reference  source  is  nearly  as  effective  in  restoring  the  imaqe 
of  an  object  at  the  height  of  the  reference  source  as  it  is  in  restorinq 
an  object  at  a  much  greater  height.  High  altitude  turbulence  is  somewhat 
more  effective  in  preventing  the  restoration  process  than  is  ground- 
level  turbulence. 


Further  work  needs  to  be  done  in  examing  the  polynomial  series 
representation  employed,  in  using  smoothing  in  the  restored  imaging 
spectrum,  and  in  extending  quoted  results  to  other  cases. 
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